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' rections to the transport coefficients of entangled polymeric systems are tested in 
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comparisons with various experimental data. It is found that quantitative descrip- 
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' tions of the viscosities, rj, dielectric relaxation time, , and diffusion coefficients, 

^ ■ D, of polymer melts can be achieved with two microscopic structural fit parameters 



whose values are in the range expected from independent theoretical or experimen- 



Q I tal information. An explanation for the (apparent) power law behaviors of r], , 



and D in (chemically distinct) melts for intermediate molecular weights as arising 
from finite size corrections, mainly the self-consistent constraint release mechanism, 
is given. The variation of tracer dielectric relaxation times from Rouse to reptation- 
like behavior upon changes of the matrix molecular weight is analyzed. Self and 
tracer diffusion constants of entangled polymer solutions can be explained by the 
theory as well, if one further parameter of the theory is adjusted. The anomalous 
scaling of the tracer diffusion coefficients in semidilute and concentrated polystyrene 
solutions, D ~ N^'^-^, is predicted to arise due to the spatial correlations of the 
entanglement constraints, termed "constraint porosity". Extensions of the theory to 
polymer tracer diffusion through polyvinylmethylether and polyacrylamide gels pro- 
vide an explanation of the observation of anomalously high molecular weight scaling 
exponents in a range where the size of the tracer, Rg, already considerably exceeds 
the gel pore size, ^g. 
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1. Introduction 



The microscopic polymer mode coupling (PMC) theory connects the dynamics of entangled 
polymeric systems to the underlying equilibrium liquid structure. This description therefore 
naturally includes corrections arising from the finite molecular weights of tracer or matrix 
polymers. In the preceding paper^, referred to as paper I, the PMC predictions with finite 
size corrections have been worked out for the transport properties of entangled solutions and 
melts, and for polymer tracer motion in gels. In the present paper these results are compared 
with data from many experiments. As the PMC approach has been detailed in the preceding 
paper^, it shall be summarized only shortly in this section and then compared to alternative 
approaches. The resulting formulae of paper I will be referenced by the equation numbers 
of that paper with a prefix I. 

The PMC approach of paper I treats the Rouse-entangled crossover problem in a simple 
interpolative manner. The PMC contributions capturing the entanglement corrections are 
combined with the Rouse, unentangled quantities, see eqns. (1.52) and (1.68) for D, eqns. 
(1.55) and (1.67) for rj, and eqns. (1.58) and (1.70) for t'^ . Appreciable deviations from 
the reptation-like asymptotes for high molecular weights are therefore not explained by 
corrections of the Rouse dynamics to the asymptotic, entangled motion. In stark contrast 
to the contour length fluctuation model of Doi^'^, the repton model of Rubinstein"^, and 
the Rouse-fluctuating-chain-end calculations of O'Connors and Ball^, the discrete nature 
of the tracer chain or its chain ends are not the cause of the deviations from the asymptotic 
molecular weight scalings. 

Finite size corrections to the transport coefficients within PMC theory result from the 
spatial and temporal correlations of the matrix constraints, called entanglements. These 
constraints are modeled with a priori reasonable assumptions. First, the local density vari- 
ations enter through a non-homogeneous compressibility, eqn. (1.45). Second, only a finite 
amphtude of the density fluctuations, or of the shear stress amphtudes ( see section 3.B of 
paper I) contributes to the elastic mesh formed by the entanglements. This amplitude is 
modeled in analogy to the entanglement plateau observed by neutron scattering from single 
polymers^^^. Except for an overall (small) amplitude, it is characterized by the entangle- 
ment length, b, eqn. (1.46). These two effects are summarily termed "constraint porosity". 



Further, the temporal decay of the matrix constraints during the disentanglement step of 
the matrix, see eqns. (1.47) and (1.48), is considered self-consistently; this process is termed 
"constraint release mechanism" and provides another parallel channel for the relaxation of 
the constraining excluded volume forces. 

Neglecting the constraint release and constraint porosity effects, PMC theory predicts 
smooth crossovers from Rouse to entangled dynamics, as can be seen from eqns. (1.52), 
(1.55), and (1.58). Rcptation-hke molecular weight scaling exponents are observed. Still, it 
is important to note that these results qualitatively differ from the corresponding reptation/ 
tube results in the prediction of non-universal asymptotic prefactors. The effects on Marko- 
vian transport coefficients of the internal mode Rouse dynamics for finite tracer molecular 
weights, which are the topic of the contour fluctuation, repton and chain-end fluctuation 
models^"^ for finite N corrections to the terminal relaxation time, are neglected within the 
PMC approach of paper 1. Recent dielectric measurements of Adachi and coworkers^' con- 
vincingly support the latter approach. A change in the molecular weight dependence of 
the end-to-end-vector relaxation time of a tracer polymer is found upon changes of the 
molecular weight of the matrix polymers^'^°; reptation-like scaling, ~ M^, is observed 
in immobile matrices of high molecular weight polymers. This effect clearly cannot be ex- 
plained by single polymer approaches but must result from considerations of the dynamics 
of the matrix constraints. 

The importance of the constraint release contributions in the viscosity has also been 
stressed by a recent simulation which focuses on the configurational rearrangements of the 
entanglement points^^. The consideration of the dynamics of the surrounding matrix is 
an aspect of PMC theory it shares with the constraint release models of Grassley^^ and 
Klein^^'^^. Whereas in these two approaches, the diffusion coefficients of the reptative and 
constraint release dynamics are phenomenologically added, in the PMC approach the effects 
of both decay channels on the entanglement friction functions depend on length scale and are 
derived from one set of approximations for the microscopic intermolecular forces, as shown in 
eqns. (1. 11) to (1.27) of paper I. Moreover, the argument of self-consistently determining the 
matrix disentanglement time from the dynamics of a single matrix polymer is not included 
in the former phenomenological constraint release formulations, whereas in PMC theory it is 



captured via eqns. (1.64) to (1.66). It appears that the plasticization model^, or the double 
reptation formahsm of Des Cloizeaux^^"^^, can be viewed as a phenomenological attempt 
to incorporate these effects. In both models, the decay of the conformational dynamics or 
stress is accelerated by taking a power of a single polymer correlation function. From the 
perspective of mode coupling approximations, non-integer powers as advocated in Ref. 5 
seem difficult to motivate microscopically. 

In des Cloizeaux's work, additional phenomenological concepts such as "time-dependent 
reptation" or "a displacement dependent diffusion constant" lead to the appearance of a 
new intermediate time scale of rather unclear physical origin. This elaboration of the tube 
model results in a stronger-than-reptation scaling of the internal relaxation time which is 
tightly connected to the anomalous frequency dependence of the high frequency wing of 
the disentanglement process in the shear modulus^^'^^. In experimental studies of shear 
moduli^^^^^ and dielectric spectra^^^^^, the high frequency wing of the disentanglement pro- 
cess consistently exhibits a power law frequency dependence, G"{ujTd ^ 1) ~ 0J~^ and 
e"{ujTjy ^ 1) ~ uj~^\ with exponents clearly much smaller than the reptation prediction, 
x,x' — 0.2 — 0.3 < 1/2. In this respect it is an important finding of the dielectric mea- 
surements that the initial decay of the disentanglement process is still more rapid, and 
consequently its high frequency wing still possesses a more shallow slope than the reptation 
prediction, even in the limit where a reptation-like mass scaling of the internal relaxation 
time is observed^'^°. Therefore, we believe a consistent and convincing physical mechanism 
for the anomalous frequency dependence of the initial decay of the disentanglement process 
remains to be identified theoretically, and it must be different from the mechanism leading 
to non-rcptation-like scaling of the longest relaxation time or viscosity. In PMC theory, 
this mechanism is the tracer shape fluctuations, as has been shown in section 3.A of paper 
I, and will further be discussed in section 3.E where the dielectric measurements of Adachi 
and coworkers^ are analyzed; see also Refs. 25 — 27. 

Whereas various different non-selfconsistent phenomenological formulations of the con- 
straint release mechanism^^^^'''^^ have been considered as extensions of the reptation/ tube 
approach, the constraint porosity effects of the PMC description, captured in e.g. the eqns. 
(1.68) and (1.69) have no analog in the phenomenological models for polymer melts and 



solutions. This follows since these models do not relate the entanglement constraints to the 
static structure of the polymeric liquid. The measurements of tracer diffusion constants in 
strongly entangled polymer solutions^^^^^ find anomalously strong molecular weight depen- 
dences, e.g. D oc N'"^-^. This rules out the constraint release mechanism as the only source 
of deviations from the reptation-like scaling. Within the reptation/ tube approach only 
an extension of the repton model, i.e. the numerical cage model of Deutsch and Madden^°, 
claims to predict molecular weight scaling of the diffusion coefficients and viscosities in agree- 
ment with the measurements in solution^^. However, from this phenomenological model it 
is not apparent why such a M-scaling is not observed in melts. Moreover, this model would 
again fail in explaining the reptation-like scaling of the internal relaxation time in the limit 
of immobile matrices, as observed by dielectric spectroscopy^'^°. A similar concern applies 
to the ideas of Rubinstein and Obukhov^^ where a close connection of the apparent power 
laws of the diffusion coefficient and viscosity, Drj oc i?^, is postulated for solutions but not 
for melts. Further non-reptation-like theoretical approaches exist^^^^^, but generally have 
not been worked out in such detail to allow comparison with the mentioned experimental 
observations and the issues of finite N corrections and apparent crossover scaling laws. 

The static correlations of the matrix constraints as described by PMC theory naturally 
generalize to amorphous media, especially gels, which are fractal on intermediate length 
scales. In the context of polymer motion through fractal media it has been recognized 
by Muthukumar and Baumgartner that strong deviations from the reptation-like behavior 
result when the medium is characterized by a length scale or mesh of the order of the size 
of the diffusing polymer^^^^^. The entropic barrier model thus agrees with PMC theory in 
predicting finite size corrections if the spatial length scale of the surrounding medium is 
not much smaller than the radius of gyration of the tracer. However, in PMC theory the 
physically reasonable but difficult to quantify concepts, such as "entropic traps" and "lack 
of topological correlations"^^ can be avoided, as the tracer dynamics is directly connected to 
the static structure of the gel and to the tracer-gel interactions. Also, the limits of small and 
infinitely large tracer sizes and smooth crossovers between regimes are naturally included 
in the PMC description but not in the entropic barrier model. In the PMC approach, 
the gel is assumed to exhibit pores of any (mesoscopic) size smaller than the gel length 



scale, ^g, which therefore determines the range of the spatial structural correlations. If 
flexible or non-rigid gels structures are considered, then another length characterizing the 
elastic stresses enters the PMC description, i.e. the eqns. (1.79) and (1.80). In contrast 
with the entropic barrier^^ considerations and its enthalpic trap generalizations^^ the PMC 
description of polymer transport through gels predicts strong finite size corrections to tracer 
diffusion even if Rg ^ ^g. 

In the present paper the PMC results of paper I shall be compared to experimental data 
in order to test the theoretical predictions. A priori estimates of the two non-universal 
parameters determining the finite size corrections to the PMC dynamics are presented in 
section 2. Section 3 introduces the experimental systems and studies examples for melt 
viscosities (sections 3.B and 3.C), melt self diffusion coefficients (section 3.C and 3.D), melt 
tracer diffusion constants (section 3.D), dielectric relaxation times in melts (section 3.E), 
transport coefficients in solutions (section 3.F), and tracer mobilities in gels (section 4). The 
paper concludes in section 5 with a summary of the results of this and the preceding paper 

1, and suggestions for future experiments to further test our predictions. 

2. Estimates of Parameters 

The PMC theory of paper I predicts that the nonasymptotic corrections to the dynamics of 
entangled polymers depend on two material and thermodynamic-state-dependent dimen- 
sionless parameters of clear microscopic meaning: S — ^p/h and a. The first parameter arises 
from porosity or correlation effects which weaken the entanglement constraints. It involves 
two measurable quantities: the polymer density-density screening length and the entangle- 
ment mesh length scale ("tube diameter" in reptation theory). Although such length scales 
can be estimated based on liquid state integral equation theory or via computer simulations, 
the best source of information is experiment. The second parameter enters the pref actor of 
the asymptotic scahng of the diffusion constant with N, i.e. D — D^{Ne/N)X^, where 
is the bare Rouse value and Xd — 32/3a. Thus, a quantifies an effective crossover degree 
of polymerization for the asymptotic behavior of f3D(o N^/N'^ where = {3a/32)Ng. 
A microscopic connection with the mean square force exerted on the probe chain center- 
of-mass by all the surrounding polymers has been derived: Xd — 32 /3a oc {\F\^)/gm, 



where Qm is the segmental number density. Hence, a is inversely proportional to the 
strength of the constraining forces (entanglements) per unit matrix number density. Based 
on liquid state theory, and the Kirkwood superposition approximation, one can derive^"^'^^ : 
{\F\'^) / Qm oc {gdipY-i where is the contact value of the interchain segment-segment radial 
distribution function which quantifies the local hard core excluded volume forces between 
chains. Although experimentally accessible in principle (and definitely via simulations), 
and gd do not appear measurable in practice. Thus, we are forced to make a priori 
estimates of a based on theoretical input. 

In this section, we discuss the available experimental information, and theoretical esti- 
mates, of the required microscopic parameters. In section 3 we treat them as adjustable 
parameters determined by fits to viscosity and/ or diffusion constant data of specific poly- 
mer solutions and melts. The independent, a priori estimates of this section provide strong 
constraints on the physically acceptable values these parameters can assume. 

We also note that, in principle, the entanglement degree of polymerization, A^g, and the 
polymer material parameters determining its bare (Rouse) dynamics, are also required as 
input to the theory. However, we always consider ratios of transport coefficients relative to 
the unentangled Rouse values, rj/ri^ and D/ D^(Nc). Thus, we employ the PMC prediction-^ 
that these ratios depend on degree of polymerization only via N/N^, where Ng is taken from 
experiment and data is plotted in the standard reduced variable format. Thus, comparison 
to experiment requires only specifying the two nonuniversal parameters 5 and a. 

A. Experimental Estimates 

Consider first the ratio 6 = ^p/b. In the Table we list rheologically-extracted (via the relation 
Gn = QraksT^cr /hY) values for the entanglement mesh, 6, of several systems in the melt"^^. 
Values appropriate for solutions containing a fraction $ of polymer follow from the known 
scaling relations^^'^^: h oc where v f=i0.75 (good solvents where oc 6 and a oc $~^/^) 
or 1/ fa 2/3 (theta solvents^^). Remarkably, these scaling relations appear to apply over a 
wide polymer concentration regime from semi dilute up to the melt^°. The density screening 
length has been measured for only a few polymers, and up to $ ~ 0.3 — 0.4 at best. For 
polystyrene, the average of many experiments yields^^: ^ 2.5A^~^ '^'^ (good solvents) and 



6.^$~^ ° (theta solvents). Of course, as the melt state is approached, theta and good 
solvents become equivalent, and hence these laws must merge into a single concentrated 
solution behavior. Elementary dense liquid considerations require that such simple scaling 
laws for the screening length cannot hold up to the melt (where < a). The screening 
length (and hence the osmotic pressure) decreases (increases) much more rapidly than in 
semidilute solution and in a non-power law fashion^^"^^. Some nonuniversality is expected 
in the concentrated and melt regime. 

Extrapolation of polystyrene semidilute good solvent scaling laws up to the melt suggests 
a screening length of roughly 3-5 A, which is expected to be an upper bound. Interestingly, 
as is seen from the Table, this extrapolation is in good accord with a value calculated based 
on the invariant packing length^^, p — The packing length rigorously emerges as 

essentially equal to C,p = {7r/3)p from the simplest analytic version of PRISM theory based 
on the idealized Gaussian thread or string models^^'"^^'^^^^^'^^. In any case, for polystyrene 
an upper bound for the parameter 5 — ^p/b ~ 0.05 in the melt can be estimated from 
the experimental data. It is probably smaller, perhaps by a factor of 2 or 3. This latter 
estimate follows from the thread PRISM based connection between the density screening 
length and the dimensionless measure of long wavelength density fluctuations^^'^^'^^"^^'^^: 
'S'o = QmkBTK = 12(^p/cr)^, where k is the isothermal compressibility and a the statistical 
segment length. From scattering experiments or thermodynamic measurements on dense 
polymer liquids, one finds that Sq is typically in the range of 0.15 — 0.3. Thus, for polystyrene 
with a statistical segment length of 7 A one obtains a ^ 1^- On the other hand, perhaps 
the dynamically-relevant mesh size for entanglements involves only the backbone mass, 
which would again result in an effective mesh size of 2 — 3 A for polystyrene. 

In solution one expects different behaviors depending on solvent quality and also chem- 
ical structure. In good solvents, where the screening length and entanglement length are 
predicted (by scaling^^ and PMC^^ theories) and measured to be proportional^^'^^, one ex- 
pects a polymer concentration-m(ie]9enden^ value of 6 but nonuniversality with respect to 
the chemical structure dependence. The latter failure of naive scaling has been emphasized 
by Adam and coworkers^^ who find a correlation of this ratio with the melt A^e- For semidi- 
lute good polystyrene solutions, from the arguments given above and known data we expect 



a larger value of 5 than in the melt, perhaps in the range of 0.2 — 0.3. Recent Monte Carlo 
simulations^''''^^ of semidilute good polymer solutions of weakly entangled chains find a value 
of 5 ~ 0.3. Even larger values are expected in theta solvents since the screening length 
increases faster than the entanglement length as polymer concentration is decreased^^'^^'^° . 
Moreover, S is expected to be polymer concentration-c?e;?enc?eni in theta solvents {S oc 
roughly) with increasing values for less concentrated polymer solutions. 

Based on the above considerations, we have a rather good idea of the magnitude of the key 
parameter 6 in the melt, and in good and theta solutions, at least for polystyrene. Although, 
we expect the concentration dependence of S is relatively universal, the magnitude may 
depend significantly on chemical structure. Direct measurements of the density screening 
length in other materials is required to determine just how large a variation is possible. 

Estimation of the parameter a, or equivalently = 32/3a, from experimental data is 
more difficult. If the asymptotic scaling regime {D oc A^~^) is achieved, then a is immediately 
determined. However, this appears possible only in melts, and not in solution. The existing 
experimental data has clearly established that Xd is significantly larger in semidilute and 
concentrated solutions of polystyrene compared with its melt value^^"^^. 

B. Theoretical Estimates 

As discussed elsewhere*^'''*^^ , theoretical estimates of the density screening length, entangle- 
ment mesh or Ne, and ga have been made based on a very simple chain model 
(Gaussian thread or string), analj^ic PRISM theory, and the RR and PMC dynamical the- 
ories. This analysis suggests or a, is nearly polymer concentration and chemical 
structure independent due to a nearly perfect inverse relationship between the contact value 
gd and the density screening length. This work provides some first principles justification for 
polymer density scalings observed experimentally, but quantitatively accurate calculations 
are difficult. Estimates based on numerical PRISM computations for solutions and melts 
which employ more chemically realistic chain models (semifiexible Koyama or RIS^^) pro- 
duce results in surprisingly good semi-quantitative agreement with the prior analytic work, 
even up to melt-like densities where g'^Sg ~ 0.1 — 0.2. Based on eqn. (1.54), we obtain the 
crude ab initio estimate oi a ^ 2/g'^So ^ 10 — 20. 



At a more detailed level, it is found that (|Fp)/^^ (proportional to a~^) decreases weakly 
with polymer concentration, and is also (modestly) quantitatively sensitive to local chemical 
features such as the effective chain aspect ratio^^'^^. Since oc {gd^pY, this behav- 

ior reflects the fact that the density and chemical structure dependences of the interchain 
segment-segment contact value and density screening length are not precisely inversely pro- 
portional to each other with a universal prefactor. This conclusion seems unavoidable to 
some degree in any concentration regime where the physical requirement for scaling ideas to 
^-PPly^ 3> CT or d, docs not hold. Thus, based on these theoretical studies, we expect some 
(modest) chemical structure variability of a (or Xd) and a decrease (increase) in magnitude 
as the melt is diluted with solvent. 

As discussed previously by comparing PMC theory with the reptation/ tube theory 
in the asymptotic limit the parameter a is easily shown to be approximately 3. Of course, 
this is only obviously relevant in polymer melts since it is the only case where the asymp- 
totic A'"-scaling of the diffusion coefficient is observed. In the reptation/ tube picture, a is 
predicted to be universal, independent of chemical structure, polymer concentration, and 
solvent quality. Such universality of a results in a predicted concentration, solvent qual- 
ity, and chemical structure independence of the constraint release corrections for viscosity 
and terminal relaxation time. Although not literally true for PMC/ PRISM theory, these 
reptation/ tube predictions of a relative insensitivity of a seem in rough accord with our mi- 
croscopic approach. We note that computer simulations have never attained the asymptotic 
regime, so they provide no knowledge of the parameter a. 

We mentioned some potential limitations of the equilibrium structural model employed in 
paper I. The simple Ornstein-Zernicke form adopted for (see eqn. (1.45)) is adequate for 
semidilute and moderately concentrated solutions. However, as the melt state is approached 
the density screening length becomes so short that it is not obvious if such a form applies. 
However, the simple equilibrium model still seems justified since experimentally the inequal- 
ities Rg,b ^ C,p apply to entangled systems"^^'^^'^^, and hence we only need to know 5*^ for 
k^p <c; 1. Another caveat is that we are treating a two component polymer-solvent mixture 
as an effective one-component system. Although this is standard practice in coarse-grained 
polymer theories, it may be inadequate for accurately computing the chemical structure. 



temperature, and concentration dependence of the structural properties required by our mi- 
croscopic dynamical theory. This fact, along with our neglect of hydrodynamic interactions, 
makes quantitatively-reliable a priori theoretical estimates of the a parameter especially 
difficult in solution. 

Finally, we point out that both a and (5 will in general depend (presumably weakly) 
on temperature via the implicit dependences of the relevant length scales, density, and 
equilibrium mean square forces. 

3. Comparison with Melt and Solution Experiments 
A. Experimental Parameters 

In order to test the description of finite size effects within the PMC theory a number of exper- 
imental data sets are studied. The shear viscosities of melts of polybutadiene^^ (PBD) and 
hydrogenated polybutadiene^^ (PBDh), and the melt self diffusion coefficients of^^ PBD, 
PBDh®^, polyisoprene^'^ (PI)i polydimcthylsiloxanc (PDMS), polyethyleneoxide®^ (PEO), 
and of polystyrene^^ (PS) are considered. In order to study the importance of the constraint 
release mechanism in melts, measurements of PS tracers in PS matrices of different molecu- 
lar weights^^, and relaxation times from dielectric spectroscopy of PI tracers in PBD melts^, 
and the neat PI melt case, are analyzed as well. Self diffusion, N — and infinite-matrix- 
molecular-weight-tracer diffusion coefficients, P ^ N, for three different concentrations of 
polystyrene dissolved in dibutyl phthalate, a good solvent, are also included^^^^^. Here and 
throughout the paper, P denotes the molecular weight or degree of polymerization of the 
matrix polymers, and M or N the tracer values, respectively. Measurements extending to 
high reduced degrees of polymerization, n = N/Ne, are of special interest. The measure- 
ments of Nemoto and coworkers on semidilute and concentrated solutions of polystyrene 
are also included. At present, no convincing theoretical understanding of these data exists 
within the reptation/ tube approach^^. 

Material parameters relevant to the entangled polymer dynamics are listed in the Table. 
The values of the molecular weights of entanglement, Mg, are the ones reported in the quoted 
experimental studies and used in our analysis. From the neutron scattering and rheological 



data survey of Ref. 46 values for the entanglement length, 6, and for the so-called packing 
length p (and estimated density screening length ^p) are included. Within Gaussian thread 
PRISM theory, is rigorously identified with the packing length, p. The data sets were 
shifted vertically on logarithmic scale in order to determine the pref actors ksT/C^Nf. or rj^ 
which lead to overlap with the Rouse result, D^oN^/kBT ^ 1/n and r]/-)]^ 1, respectively, 
for small molecular weights, n = N/Ng = M/M^ <C 1. Note that in all cases data for n < 1 
were available, and as little arbitrariness was introduced in the shifts as possible. The PMC 
dynamical parameters, a and 6, listed in the Table follow from the fits discussed in the 
following sections, and carry significant uncertainties. 

B. Melt Viscosity 

A PMC analysis of the reduced viscosity, rj/r]^, versus reduced molecular weight, n = N/Ne, 
highlights the effects of the self-consistent constraint release mechanism and involves one fit 
parameter only, the inverse strength parameter a. Comparing the asymptotic reptation and 
PMC predictions for the viscosity at n — oo, the universal value a ~ 3.2 is estimated. Figure 
|l] shows the viscosities of PBD melts for which the highest experimental reduced molecular 
weights have been reported^^. The theoretical result, eqn. (1.74), has been used. The 
PMC result (at least) semi quantitatively accounts for the measured data for all molecular 
weights. The asymptotic scaling, t] ~ A^^ (not shown in figure |l]), is not yet exhibited by the 
data nor by the PMC fit even though values n = 10^ are achieved. The matrix constraints 
decay on the same time scale which determines the viscosity, i.e the disentanglement time, 

, and less friction therefore arises than in the asymptotic limit of frozen entanglement 
constraints. The constraint porosity effects, i.e. the spatial correlations of the entanglements 
characterized by nonzero values of the density screening length, ^p, and the entanglement 
length, b, have been neglected in figure as they do not noticeably affect the shape of the 
theoretical curve for expected values of the length scale ratio in the melt, 6 = ^p/b ~ 1/20. 

A second fit of the PMC results to the data is shown in figure 0. It assumes, as has been 
pointed out by O'Connor and Ball^, that the monomeric friction coefficient, (q, found in 
the Rouse regime in Ref. 61 is somewhat too high. O'Connor and Ball argue by comparing 
the friction coefficients from Ref. 61 with an earlier measurement by Roovers^^, that the 



Rouse results of Ref. 61 are a factor 2.346 too high, Co — 2.346Co- If this correction of 
the monomeric friction coefficient is taken into account, by fitting the expression T]/ri^ = 
1 + 2^64'^'^^ data instead of eqn. (1.74), then a quantitative fit of the PMC resuhs to 

the measured viscosities becomes possible in the full range of molecular weights excluding 
the highest data point only. The value found for the strength parameter a, a ~ 7.7 without, 
and a ~ 4.9 with, this correction, then lies closer to the reptation prediction of a ~ 3 and 
also to values found in other polymer melts, e.g. PBDh (see below). Such values of a are 
also expected from the observation of the apparent power rj ~ A^^'^, as demonstrated in the 
model calculations shown in section 5. A of paper I, where a power rj ~ A^^'^ is observed for 
values of a below 5. 

After correction of the monomeric friction coefficient, all data except for the highest 
molecular weight point (which is not reliable^^) fall close to the PMC fit. This is especially 
apparent from the inset in figure ^ where the relative variation of the viscosity compared to 
its asymptotic behavior, 7]/M^, is shown. The data scatter rather unsystematically around 
the fit with deviations well compatible with the error bars reported in Ref. 61 due to 
uncertainties in the measurements of the molecular weight and the viscosity. The slow 
approach of the constraint release corrections to the asymptote, eqn. (1.69), explains why 
the asymptotic rj/M^ const, is not yet observed. 

From figure and especially from the inset, it is suggested that the apparent power law 
behavior of the viscosity in an intermediate molecular weight range, rj ~ A^^'^, is not the sign 
of a true power law scaling behavior but arises from the competition of a finite size effect, 
namely the self-consistent constraint release mechanism, and the asymptotic, reptation-like 
scaling, t] ~ M^. Clear observations of the true asymptotic power are not expected in the 
experimentally accessible molecular weight range, as the constraint release correction factor, 
P{n/a'^), in eqn. (1.67), shows an extremely slow approach to unity, 1 — /3 — > 4/3(a^/n)^/^ 
eqn. (1.73). Effectively, the slow variation of (3 around n ^ 10^ may mimic the asymptote, 
eqn. (1.54), which, however, correctly describes the data for much higher rescaled molecular 
weights n only. 

Unfortunately, no diffusion data exist for PBD melts out to such high molecular weights 
as used for the viscosity measurements. In the next section, the combined comparison of 



the theory with diffusion and viscosity data provides a further more stringent test of the 
physical picture of PMC theory. 

C. Combined Viscosity and Self Diffusion Melt Data 

Of special value to test the close connection of the center-of-mass and conformational dy- 
namics as it is predicted by PMC theory are data on both dynamics for identical polymer 
systems. An interesting example are the measurements in PBDh by Pearson and coworkers 
which extend to rather high degrees of polymerization and provide self diffusion constants 
and viscosities for the identical samples^^. The viscosity analysis, using the PMC result, 
eqns. (1-74), fixes the inverse strength parameter a and therefore the constraint release 
contribution to r/ and D^. The value a = 5.2 gives the best least square fit to the viscosity 
data and is rather close to the value fitted to the PBD viscosities. Figure || displays the data 
and PMC results. A universal value of a is not expected by PMC theory; the magnitude 
agrees well with the PRISM estimate and the universal value predicted by reptation theory. 
The apparent power law behavior of the viscosity^^, rj ~ A^^'^, is again explained as a finite 
size crossover effect due to the constraint release mechanism. 

A slightly smaller parameter, a = 4.3, provides the combined best PMC fit to the 
viscosity and the self diffusion coefficient as a function of reduced molecular weight. The 
lower part of figure ^ compares the two viscosity fits which exhibit similar fit quality. In the 
upper part of figure |^ the self diffusion coefficients are seen to overshoot slightly compared 
to the reptation like asymptote, D D^/Xdu, eqn. (1.52). This small enhancement of the 
diffusion constant again results from the constraint release mechanism as is shown by the 
PMC fit using eqns. (1.64,1.68,1.69) and the same parameter, a = 4.3, as the viscosity fit. 
The small density screening length, 6 = C,p/b ~ 0.03, found in this analysis indicates small 
constraint porosity contributions to D. Moreover, this value for 6 is in good accord with the 
independently measured length scales listed in the Table. An independent least square fit 
to the diffusivities alone leads to a slightly smaller value of a = 3.5, but somewhat a larger 
constraint porosity contribution, 6 = 0.04. The rather small variations of the parameters, 
a = 4.3 ± 1 and 6 = 0.03 ± 0.01 for PBDh, however, are well compatible with both the 
simplifications of the polymer matrix model, eqns. (1.45) and (1.46), used in the theoretical 



approach, and a priori estimates. 

It is the decay of the matrix constraints which is the dominant finite size eff'ect in poly- 
meric melts. Therefore, the larger finite size corrections in the viscosity compared to the 
diffusion coefficients are expected in PMC theory. In paper I, it was discussed that the 
larger intramolecular weighting of the matrix entanglement contributions for long wave- 
lengths leads to larger constraint release effects in the conformational dynamics than in 
the center-of-mass motion. The observation in polymer melts of the asymptotic scaling, 
D ~ A^~^, in the self and tracer (for P > N) diffusion coefficients, but not in the shear 
viscosities, is therefore naturally explained by the different intramolecular correlations cap- 
tured in the S and M memory functions of PMC theory. See sections 4 and 5, and figures 
4 to 8, in paper I for more discussions of this connection. 

D. Melt Tracer and Self Diffusion Constants 

The small corrections of the finite size effects in the melt diffusion data have two origins 
in PMC theory. First, the constraint release mechanism is rather ineffective for the center- 
of-mass motion because of the local nature of the S memory function. Second, in melts 
the spatial correlations of the matrix constraints are characterized by rather small lengths 
scales, and b, and a length scale ratio 6 = C,p/b, much smaller than in the solution cases. 

Self diffusion data for different melts do not fall on a single master curve, D/D^ = 
f{N/Ne), as predicted by reptation, because of the variations in the (inverse) entanglement 
strength factor, a. Independent fits to a number of polymer melt self diffusion constants 
(PBD, PBDh, PDMS, PEG, PI, PS), however, show that values of a fall in a rather tight 
range of a = 2.5 — 5. The density screening lengths, ^p, and length scale ratios, 6 = C,p/b, 
found in these fits lie consistently below 6 = 0.06 as expected for dense melts (see the Table). 
Figure ^ displays a representative set of self diffusion data and corresponding PMC fits. 

A particularly interesting set of self diffusion data is included in figure |^. Diffusion 
constants of PDMS^^ exhibit an apparent power law, D ~ A^~^-^, up to ca. log{N / N^) = 1.2. 
Such high apparent power law exponents for rather small degrees of polymerization are well 
rationalized by the PMC finite size effects, mainly the constraint release mechanism, as 
shown by the PMC fit in figure ^. During this fit three parameters, a, Xd and 6, were varied 



freely in eqns. (1.64,1.68,1.69). Supporting the close connection of the center-of-mass and 
the conformational motion predicted by PMC theory and leading to eqn. (1.53) for Ad, the 
two resulting fit parameters, a and A^), almost satisfy eqn. (1.53), 32/3a ~ 3.8 compared to 
\d = 3.6. 

A common PMC fit to a number of self diffusion data sets using eqns. (1.64,1.68,1.69) is 
shown in figure ^. Rather small enhancements of the self-diffusion constants compared to 
the asymptotic value result from the small consequences of the constraint release mechanism 
in the melt. The small value of the extracted 6,6 = 0.05, indicates that constraint porosity 
can be neglected in this fit. 

Without showing a corresponding figure, let us report that for tracer polymers in a 
strongly entangled polymer melt, where the matrix molecular weight far exceeds the tracer 
molecular weight, P ^ N, almost no deviations from the asymptotic prediction, D = 
D^/{1 + Xon), eqn. (1.52), are observed in the PMC results in agreement with experimental 
measurements^^. In figure 7 of paper I, model calculations for melt-like parameters, a = 3 
and 6 < 0.05, show a smooth crossover of this infinite matrix molecular weight tracer 
diffusion constant from its Rouse, D ~ A^~^, to asymptotic, D ~ A^~^, behavior with no 
overshooting or higher apparent exponents in intermediate N regimes. This behavior clearly 
contrasts with the rather strong overshooting or rather high apparent exponent, D ~ A^~^'^, 
caused by the constraint release mechanism in the PDMS self diffusion data of figure |^. 
Note that for n = N/Ne = 4.4, the PDMS self diffusion constant lies a factor 2.1 above the 
asymptote obtained by neglecting finites size effects. The individual analysis of the PBD 
and PEO self diffusion coefficients in figure ^ further support these findings; small length 
scale ratios are observed and dominating constraint release corrections to the asymptotic 
behavior. 

The importance of the dynamics of the matrix constraints on the tracer diffusion coeffi- 
cients can most comprehensively be studied from a set of data for different tracer, M, and 
matrix polymer, P, molecular weights. Recalling that our approach applies to entangled 
polymer matrices only, due to the approximations in the RR model, we study the data of 
Green and Kramer^^ for PS tracer in PS melts for P > P^. As the authors did not publish 
tables of their extended set of data we use the reported fitted theoretical curves in order to 



compare with the PMC results. Note that Green and Kramer could fit their data equally 
well with the models of Grassley^^ and Klein^^. Both models mainly differ in the steep- 
ness of the apparent P dependence of the tracer diffusion constants in the crossover range, 
P < N. We use the fitted model results of Grassley because we expect easier agreement of 
PMC theory with the model by Klein; in figures 8 and 9 of paper I, P-dependences around 
Klein's prediction, D ~ P~^'^, are seen in the PMC results. Grassley 's expression exhibits 
a stronger P-dependence, 

D(M, P) « ^(1 + a^^MlMP-^) . (1) 

Note that eqn. (1) does not reproduce the Rouse model for small tracer molecular weights, 
M < Me. Therefore we normalize it to unity for M = M^, i.e. study D{M, P) / D{Me, P), 
and do not require the connection of a and Xd, eqn. (1.53), i.e. the relative ratio of Rouse 
to asymptotic PMC result, in the PMC fit. The experimental parameters, Mg as in the 
Table and ttcr = ||^(^)^~^ with z = 3.5, are taken from Ref. 66 and therefore mimic 
the experimental data for the PS tracer in PS melt systems. The values for the tracer and 
matrix polymer molecular weights used in figure ^ span the experimentally studied range. 
Even though the PMC results in figure ^ do not reproduce the D ~ P^^ prediction of 
Grassley for the tracer diffusion constant in the the constraint release dominated region, 
overall a reasonable fit to Grassley's model can be achieved as shown. The value a = 3 
is chosen a priori as appropriate for polymer melts, and the parameters 6 = 0.07 and 
Xd = 1-58 are found from least square fits. Comparing the shapes of the PMC results with 
the actual diffusion data of Ref. 66 qualitative agreement can be deduced. It appears that the 
PMC theory adequately describes the constraint release decay of the matrix entanglement 
constraints as observed in tracer diffusion measurements in entangled polymeric melts. 

E. Dielectric Relaxation Times 

An especially powerful technique to study the conformational tracer dynamics, especially 
the end-to-end vector fluctuations, is the dielectric spectroscopy of type-A tracer polymers 
in type-B polymer melts^'^°. Type-A polymers are defined to possess monomeric dipole 
moments parallel to the chain contour, whereas type-B polymers possess perpendicular 



moments only^^. Adachi and coworkers have studied such a system, PI tracer in PBD 
melts, for varying tracer and matrix molecular weights^. They report the molecular weight 
dependence of the tracer dielectric relaxation time and the high frequency asymptote of the 
disentanglement process in the dielectric spectrum. It is important to recall from paper I that 
PMC theory identifies two different mechanisms, one which is responsible for the initial decay 
of the end-to-end vector correlation function, {P{t) ■ P(0))/(P(0) • P(0)) - 1 ~ (t/NHo)'', 
and two other ones, which affect the dielectric relaxation time. At first, let us discuss the 
latter, namely the constraint release and constraint porosity effects, and compare them with 
the experimental data for , the dielectric relaxation time. 

As shown in section 5.D of paper I, PMC theory predicts three different asymptotic 
behaviors for the relaxation time of the end-to-end vector correlation function of a tracer 
polymer. In unentangled matrices, p = P/Pe < 1, the Rouse behavior, ~ A^^ should 
be observed. In strongly entangled matrices, where the surrounding polymers are immobile 
relative to the tracer, the dielectric follows the internal disentanglement time, oc r^, 
and exhibits reptation-like scaling, ~ A^^. For tracers which are much larger than the 
matrix polymers, which are entangled however, i.e. n — N/N^ » p » 1, PMC theory 
finds an intermediate scaling behavior, ~ N^p^^^n^^^^ . Note that these results, more 
strongly than our results for the internal or the center-of-mass dynamics, depend on the 
approximations made when diagonalizing the equations of motion, eqn. (I.l), where chain 
end effects have been neglected^^. 

Also, for a PI tracer in a PBD matrix system-specific polymer interactions could 
strongly influence these results, especially the parameters characterizing the tracer dynam- 
ics. Whereas the parameters, Pe, a"*, and 6, specifying the dynamics of the entangled melt 
correspond to the ones found (for example) from independent viscosity measurements, the 
tracer parameters, N^. and a*"^, must differ from the corresponding values found for the one- 
component system of tracer polymers. As eqns. (1.35) and (1.54) show, these parameters are 
determined by the tracer-matrix interactions. Moreover, the importance of the constraint 
release mechanism depends on the ratio of the monomcric time scales of the tracer and 
matrix polymers. The corresponding variable is r = r^^Pe/r^Ne in eqn. (1.64), where Tq is 
the monomeric time scale of the tracer in the specific matrix it is embedded in, not in its 



own melt state. 

Adachi and coworkers have observed that the dielectric relaxation times, , of PI 
tracers in PBD matrices measured in the molecular weight ranges, —0.2 < log^g^ — 1-45 
and 0.0 < log^gP — 2.03, can be fitted with the following formula: 

logio(r7iVe'ro) = -1.92 + (3 - e-(^/«^°°)"') logio(M/1900) , (2) 

where tqN^ = 10^^'^^sec, and the molecular weights of entanglement are Pe = 1850 and 
Me = 5000. Note that the range 0.03 < log^^Q^ < 1.03 is of special interest, as there the 
effective exponent is found to lie between the expected asymptotic behaviors, ~ A^^ or 
~ A^'^. These limits. Rouse for P < P^, and rept at ion-like for P ^ Pg, are also predicted 
by PMC theory. For intermediate matrix molecular weights, the PMC approach does not 
predict simple power laws as reported^ in eqn. (2). A crossover from ~ A^'^ for N <^ P, 
to ~ ]V^^/^ for N ^ P, for entangled matrices, P > Pe, results in some curvature. The 
PMC curves matched to eqn. (2) therefore do not fit quantitatively. In the comparison of 
figure P the PBD melt is a priori characterized by parameters, a = 4 and 6 = 0.05, found 
previously in sections 3. A and 3.B. Then three parameters are varied in the fitting procedure. 
For simplicity, Mg of the PI tracer in the PBD melt is chosen to be the one for a PI melt, 
and a*^ = 12 is matched to the asymptotic behavior, t'^/N^tq ~ f|^- The fitted small 
value of the time scale ratio, r = Pe/r^Nf. = 0.005, indicates that the tracer polymer is 
relatively mobile compared to the matrix polymers. As is evident from eqn. (1.64), the effect 
of a small parameter r can also be interpreted as a change in N^, O oc rN~^^'^, and without 
further independent experimental information a closer identification of the PMC model 
parameters is impossible. Also, the errors introduced into the PMC results by the neglect 
of chain end corrections^^ are unknown at present. In figure ^ the reported change^'^° from 
Rouse to rept at ion-like molecular weight scaling of is qualitatively correctly reproduced 
by the PMC results, but with quantitative deviations. Instead of a scaling exponent which 
continuously depends on matrix molecular weight^, three exponents, 2, 19/8 and 3, are 
predicted, whose observation depends on the ratios N/P and P/Pe- 

Let us recall the predicted behavior of the melt dielectric relaxation time shown in figure 
10 of paper I. There, the dielectric relaxation time of one of the matrix polymers, i.e. N = P 
and identical chemistry, was observed to scale like ~ A^^'^ in an intermediate molecular 



weight range, 0.5 < logiQ^N/Ne) <2. Figure |^ shows the dielectric relaxation times for 
PI homopolymer melts^^^^*^ and PMC results for melt like parameters a = 2.2 — 3.5 and 
6 = 0.05. Both values agree with the estimates given in section 2 and are compatible with 
the analysis of the self diffusion data^^ shown in figure ^ Thus, PMC theory also achieves 
a quantitative description of the PI melt dielectric relaxation times mainly by adjusting the 
self-consistent constraint release mechanism efffects via the paramter a. 

The molecular weight dependence of the dielectric relaxation time thus is influenced 
by the finite size size effects which result from the spatial and temporal correlations of the 
matrix constraints. The initial decay of the disentanglement process in the end-to-end vector 
correlation function, however, is determined by the tracer shape fluctuations in PMC theory. 
As shown in section 3. A of paper I, the reptation like behavior, e"{ujT^ ^ 1) ~ tu^^/^, emerges 
from the PMC description only if the internal contributions to the collective tracer dynamics 
in the RR model are neglected, which is not justifiable. Even in the case of tracer motion in 
matrices of immobile polymers, where the rept at ion-like scaling, ~ A^^, is predicted, the 
tracer shape fluctuations of PMC theory still speed up the initial disentanglement process, 
resulting in e"{uJT^ ^ 1) ~ cu^^, with 9/32 < x < 3/8 (see section 3. A of paper I). These 
predictions qualitatively agree with the observation^ of (i) frequency power law exponents 
around x = 0.22 to 0.31 for PI tracers in strongly entangled PBD matrices (A^ ^ P and 
P ^ Pe)j where ~ A^^, (ii) exponents x ~ 0.53 in the Rouse limit, P^ N <^ P, 
and {Hi) x values in between for matrix molecular weights in between. Only numerical 
calculations'''''''^ of the full PMC equations can determine whether quantitative agreement 
results. In this respect, it is important that the constraint release correction to the initial 
decay of the disentanglement process can be neglected, as it sets in only for t > r^,. 

F. Solution Transport Coefficients 

It has been rather well established that the molecular weight dependence of the viscosity 
of entangled polymeric solutions is very similar, and perhaps identical to within experi- 
mental error, to the behavior in polymer melts^°'^^. For example, in solutions of PS in a 
good solvent, Nemoto and coworkers^^'^^ found that for the semidilute and concentrated 
solutions the viscosity shows an apparent power law behavior, rj ~ jys.sio.i ^j-^g range 



< \og{N/Ne) < 1.7. As figure || shows, PMC theory can quantitatively account for these 
behaviors, and (inverse) entanglement strength factors in the physically expected range of 
a = 2.5 for 40 wt% PS and a = 4.0 for 13 wt% PS are found. The similarity of the entangled 
solution viscosity to the melt viscosities finds a clear physical origin within PMC theory in 
the similarity of a for solutions and melts. From the PRISM estimates, this density and con- 
centration independence of a has been anticipated. For the semidilute concentration data, 
13 wt % in figure |^, it is not clear how to unambiguously compare theory and experimental 
data. For low molecular weights, N < Ng, the data exhibit a slightly stronger dependence 
on molecular weight than explained by the Rouse model, which is assumed in the PMC 
approach to apply for < Ne- The vertical shift of the 13 wt % data therefore is unclear. 
We estimate that this is the origin of the variation of a between 2.5 and 4 found in the fits 
of figure these fits result if the experimental viscosities are normalized to ry/M = 1 for 

M = Me. 

The PMC theory together with the PRISM estimates of the dynamical parameters ex- 
plain that the viscosity versus molecular weight curves for solutions and melts are similar. 
Note that the theoretical predictions for the concentration dependence of N^. can in principle 
be used, although in our case Ne was taken from experiment^^^^^. The similarity of solution 
and melt viscosities results in PMC theory from the equal influence of the matrix constraint 
decay relative to the probe diffusive motion in the conformational memory function, M{t). 
Important for this result are two physical aspects. First, rather long ranged intramolecular 
correlations dominate the conformational friction function; see the discussion of J*^, eqn. 
(1.63). This leads to insensitivity of r] to local length scales like the density screening length. 
Second, the (near) cancellation of the concentration dependence in the PMC entanglement 
strength factor^^, 1/a ~ g'^So ~ g'^. Although at lower solution densities the polymer den- 
sity fluctuations, or the osmotic compressibihty, are strongly enhanced (larger 5*0), the local 
contact probabilities of monomers on different polymers, g^, diminish, thereby resulting in 
a ~ constant with varying polymer concentration. 

When studying the tracer and self diffusion coefficients of the same PS in good solvent 
systems^^^^*^, it is the constraint release mechanism which leads to the differences between 
the self, N = P, and (infinite matrix molecular weight) tracer, P ^ N, diffusion coefficients. 



In figure ^ one observes a similar enhancement of tlie solution self diffusion constants rela- 
tive to the solution tracer diffusivities as has been observed in figure § for the melt diffusion 
data. There the melt self diffusion constant at N/N^ = was a factor 2.1 larger than the 
asymptotic expression, eqn. (1.52), which neglects finite size and especially constraint release 
effects. Due to the similar constraint release parameter, a = 2.8, a priori chosen in the analy- 
sis of figure ^, a similar enhancement of the self diffusion data (affected by constraint release) 
to the P ^ A^-tracer diffusion constants (unaffected by constraint release) is observed. For 
example, at n = N/N^ = 4.4 for the 40 wt% PS solution, ^(A^ = P)/D{P > iV) ^ 2.1 is 
found in figure 

Even though the constraint release mechanism cannot be of relevance to the (infinite 
matrix molecular weight) tracer diffusion constants, they do not exhibit the reptation like 
scaling, D ~ A^~^, in dense PS solutions^^^^^. The experimental data of Nemoto and 
coworkers^^ which extend to rather high molecular weights, n < 70, exhibit an apparent 
power law behavior, D ~ j\r-2.5±o.i^ ^^^^^ jjg ^^jj below the corresponding melt data at high 
molecular weights plotted in an isofriction comparison for fixed n = N/N^.. In figure ^ the 
measured diffusion data are compared to this power law, to a PMC fit, and to the asymptotic 
PMC prediction, eqn. (1.52) with \d = 19.2. The PMC fit uses an independently, from the 
viscosity analysis, determined a, and freely varies 5 and Xr,. Note that the short dashed curve 
in figure ^ shows the asymptotic PMC result including the Rouse, unentangled behavior, 
i.e. D = D^/ (1 + Xon). The enhancement of the diffusion, especially of the tracer diffusion 
coefficients, relative to the asymptotic behavior, D ~ A^~^, therefore cannot result from the 
crossover to the Rouse law, D ~ 1/A^. The enhancement of the PS diffusion coefficients 
in entangled, PS solutions relative to the short dashed asymptote in figure § is a finite size 
effect according to the PMC fits included in figure |. 

From the true PMC asymptote, D ~ D^/Xdu, one can deduce that the asymptotic 
prefactor, Xo, which, in the comparison with the melt diffusion constants in figure was 
found to be around 3 in melts, has to be appreciably larger in solutions. A value of Ad 
smaller than Ad = 12 appears not compatible with the solution data as they exhibit an 
increased molecular weight dependence out to the highest A^ points. Such a large value 
of Xd is not compatible with the relation A^ = 32/3a, eqn. (1.53), as the value of a lies 



in between 2 and 5 as has been discussed in context with the solution viscosity and the 
difference of self to tracer solution diffusion constants. At present it is not clear how this 
violation of eqn. (1.53) can arise in PMC theory. In the discussion of section 5 we will list 
possible mechanism. One possible origin is the simplified treatment of the two-component 
system, i.e. polymer solvent mixture, as an effective one-component system. Lacking a 
quantitative theoretical understanding of A/? in solutions we suggest using it as a free fit 
parameter in eqns. (1.64,1.68,1.69), while a = 3 ± 1 may be fixed a priori from independent 
viscosity measurements. 

The PMC fits to the self and tracer diffusion constants in entangled, concentrated poly- 
mer solutions in figure ^ thus show that again the apparent power law behavior, D ~ A^~^'^, 
can be accounted for by a finite size correction of PMC theory, the constraint porosity. It 
arises from the rather large density screening length, 5 = ^p/h = 0.3, found in the fits. Note 
that the constraint porosity corrections, which cause the strong deviations of the solution 
diffusivities from the asymptote, D ~ A^~^, in figure ^ have no analog in the reptation/ tube 
picture. Consequently, no consistent theoretical understanding within that approach has up 
to now been presented for the solution data of Nemoto and coworkers^^ which explains the 
difference^^'^^ of solution and melt diffusivities and the melt results as well. 

Within the PMC theory, the enhanced effective exponent results from the gradual build 
up of the full constraint amplitudes a tracer polymer feels as it compresses the surrounding 
entanglement mesh. Only when the size of the polymer tracer far exceeds the entanglement 
length, b, and the compressibility or density screening length scale, ^p, are the full matrix 
constraint amplitudes effective and lead to the reptation like asymptotic behavior. As the 
reptation/ tube approach neglects the spatial variation of the compressibility by simplifying 
the spatial correlations of the entanglement constraints to arise solely from the entanglement 
length {6 = effectively), it misses this effect and cannot explain the solution diffusion 
coefficients. The PMC approach naturally includes the constraint porosity effects as it 
starts from the microscopic correlations of the polymeric liquid structure. This is especially 
necessary in polymer solutions where the density screening length is not small compared 
to the entanglement length, 6 = C,p/b = 0.3 is found for the PS solution data of Nemoto 
et.al.^^^^^. This length scale ratio is larger in solutions than in melts as has been argued 



in section 2. It is rather intriguing, but not understood at present, that the ratios of the 
asymptotic prefactor. Ad, in solution to the one in melts, and of the length scale ratio 6 in 
solution to its value in melts, are rather similar: A^''/Ag'^'* ~ ^soy^meit ^ ^ — q 

The possibility of a common fit with one set of parameters of self diffusion relative to 
its tracer diffusion analog at three different polymer concentrations in good solvents shows 
that the constraint release enhancement of self to tracer diffusion constants is correctly 
captured in PMC theory, and that the PRISM prediction of concentration independence of 
the parameters a and 6 in solutions is accurate. Note that variations of a around a = 3 ± 1 
do not worsen the fits. This point is addressed in figure |1^ by independent least square 
fits to the individual data sets; a = 2.83 is apriori fixed. Variations of 11 < Ad < 39 
and 0.2 < 6 < 0.5 again emphasize that both parameters are larger for solutions than for 
melts, but do not exhibit systematic polymer concentration dependences in the semidilute 
and concentrated solution regimes. 

The enhancement of the diffusion coefficients compared to the asymptotic, reptation like 
results, D ~ A^~^, is considerable even for large molecular weights. For n = N/N^ = 4.4 and 
the 40 wt% PS solution, the tracer diffusion coefficient lies a factor 6, and the self diffusion 
constant a factor 12 above their common asymptote in figure ^ at n = 30 for the 18 wt% 
solution the enhancement factors are still 2.2 and 1.9, respectively. Extending the PMC fits, 
a merging of the diffusion data with the asymptote can be expected only for even larger 
molecular weights, ca. n ^ 300. 

4. Polymer Tracer Motion in Gels: Comparison with Experiment 

As the extension of the PMC theory to polymer tracer diffusion through gels has been pre- 
sented in paper I, and since tracer motion through gels is the original problem studied by 
de Gennes^^, two experimental studies of this problem are included in this paper. The low 
field mobility, /i = DN, from electrophoresis measurements of polystyrenesulfonate (PSS) 
in polyacrylamide (PAA) gels^^ and the tracer diffusion measurements of PS in polyvinyl- 
methylether (PVME) gels'^^, are considered. 

The PMC theory can also describe the motion of polymer tracers through random ma- 
terials, where especially gels are of interest (see paper I). The physical aspects of constraint 



porosity in the PMC description naturally generalize to fractal media like gels. Two dif- 
ferent gel models were mentioned: "hard gels", eqn. (1.76), which are characterized by 
rigid density structures which pose the constraints on the tracer dynamics, and "soft gels" , 
eqn. (1.77), where the flexibility of the matrix structure results in much smaller elastic con- 
straints. Simple extensions of the PMC results for polymer melts and solutions just replace 
the tracer and matrix equilibrium structure by a fractal tracer or gel structure respectively. 
Due to the time independence of the arrested gel structure, the constraint release mechanism 
is absent in gels or amorphous solids in general. In PMC theory strong enhancements of 
the tracer diffusivities relative to the asymptotic power laws are predicted from the fractal 
generalizations of the constraint porosity effects. 

As the tracer-gel chemical interactions may show strong variations, the possibility to 
make a priori PRISM estimates for the dynamical parameters is missing at the present state 
of the art of PRISM^^'^^. From physical considerations, as the entanglement asymptotic 
prefactors, i.e. the A's in eqns. (1.67,1.68,1.79), clearly are proportional to the tracer matrix 
interactions, a^^ oc A oc (|-FtgP)/f?m where the ¥tg are the tracer-gel interaction forces, and 
as stronger chemical interactions, like preferential adhesion, clustering or partial wetting, are 
expected for polymers in gel pores, it can safely be concluded that the asymptotic prefactor 
\g in eqn. (1.79) may be significantly larger than estimated from simple homopolymer 
solution or melt systems. 

As has been shown in section 6 of paper I, the most important variable determining 
the tracer diffusivities is the ratio of gel pore size, ^g, to tracer size, Rg. Large finite size 
corrections are predicted by eqns. (1.79) and (1.80) for tracers which persist even if Rg ^ E,g 
but are not yet in the true asymptotic limit, which may require Rg/^g = 10^. As the 
asymptotic regime is reached when Rg ^ C,g, changes in the mass-size scaling of the tracer 
polymer do affect the shape of the pre-asymptotic apparent scaling region. Variations in 
the gel fractal dimension, however, mainly shift the fitted (unknown) material parameters, 
but do not significantly influence the shape of the D{N) curves. 

Figure |1^ reproduces the low field mobility data''^, fi = DN, of PSS in PAA gels of 
different gel concentrations, c. From the data in figure |lT] it is apparent that much stronger 
molecular weight dependences of /i are observed than predicted by either Rouse or reptation 



theory. This has prompted Muthukumar and Baumgartner to develop the idea of entropic 
trapping^^^^^ as an alternative transport mechanism which may be dominant in the Rg ~ 
regime.. We will analyze this observation using the finite size effect, namely the constraint 
porosity, captured by PMC theory for polymer melts and solutions. 

As the PAA gel structure is assumed to be similar to the solution structure of PAA'''^, 
the fractal dimension appropriate for polymers in good solvents. Dp = 5/3 in Flory approx- 
imation, is chosen in eqn. (1.75). Also the soft gel model clearly applies as the crosslinked 
polymer matrix does not exhibit significantly stronger constraints than the entanglement 
mesh of the corresponding solution, but only time independent ones. We use this consid- 
eration also to a priori fix the ratio of entanglement length to gel pore size to the value, 
^ = ^/b = 0.3, found for solutions in the previous section. Following the suggestion in Ref. 
74 the mass size scaling of the tracer polymer is also described by the Flory behavior, i.e. 
z/ = 3/5 in eqn. (1.74). 

Figure |ll] includes fits with the PMC results, eqns. (1.79) and (1.80), for the model of 
tracer and gel equilibrium structures described above. Least square optimization was used 
in order to determine the tracer-gel interaction parameter, A^, and the gel pore size, ^g, for 
three different gel concentrations. It appears physically very reasonable that the tracer-gel 
interactions are gel density independent, and = 54 ± 5 is found for all three densities. 
Gel pore sizes are reported in figure [TT] by stating the PSS degree of polymerization, A^, 
giving R^^^{N) = C,g. The expected decrease of the gel pore size with increasing gel density 
is clearly observed. A least square fit to the three data points yields ^g ~ c^^'^, which 
overestimates the expected exponent^^ (^g ~ c^^^'^ for semidilute good solutions) and lies 
closer to the theta solvent values {^g ~ c"^ for semidilute theta solutions). 

In contrast to the expectations of the entropic barrier model of Muthukumar and 
Baumgartner^^, the steepest region of the mobility versus molecular weight curve is found 
to lie at Rg/C,g ^ 1- The entropic barrier model leads to strongly tracer molecular weight 
dependent diffusivities when the tracer random coil just fits into the largest gel pores. Al- 
though the PMC constraint porosity effects set in when Rg ^ ^g, it extends out to very 
large tracers and leads to high apparent exponents even when Rg ^ C,g- Also in contrast to 
the entropic barrier model, the PMC approach includes the smooth crossovers to the non- 



interacting Rouse and the asymptotic generalized reptative limits. It not only describes 
the region of high effective power laws, but also describes a smooth variation for all tracer 
molecular weights. Note that, for the fractal dimensionalities chosen for PSS tracers in PAA 
gels in figure |TT], the tracer diffusion constants do not crossover to the reptation behav- 
ior, /i = ND ~ 1/iV, at extremely large tracer molecular weights, but instead exhibit the 
much weaker dependence, jj, ~ A^~^/^, as follows from eqn. (1.79). This weak asymptotic 
molecular weight dependence is compatible with the data as the crossover to this asymp- 
tote is predicted to take place at much higher (unobservable) only. In the intermediate 
range effective exponents up to /i ~ DN ~ A^~^'^ have been reported^"^, which lie in the 
range and below the upper bound deduced in paper I for these fractal dimensionalities, 
/i ~ j\j-2{2-u{d-Dp)) ^ jY-12/5^ Deviations of the fits in figure |ll] for small molecular weights 



presumably arise from the inadequacy of the Rouse model for short polymers in gels. 

Another study of the tracer diffusion of PS through PVME gels''^ was aimed at studying 
the problem originally considered by de Gennes when introducing the reptation concept 
The gels and their flexible, neutral polymer solution precursors were well characterized 
and considerably crosslinked, and should therefore be a model system for testing reptation 
ideas. It was suggested that the gels fall in the light strangulation regime, where the polymer 
strands between chemical crosslinks are somewhat shorter than the ones between the time 
dependent entanglements of the non-crosslinked polymer solution. Following the estimates 
presented in Ref. 75 we assume that for the three different crosslinking densities studied the 
length scale ratios, ^g/b are somewhat larger than 6 found for solutions, ^g/b = 0.3, 0.4 and 
0.5 are chosen as reasonable estimates. The smallest value, found for PS solutions, would 
correspond to the border of the strangulation regime, and ^g/b = 0.5 to a molecular weight 
in between crosslinks a factor of about 2.8 smaller. 

It came as a surprising finding, which was claimed to be in disagreement with reptation 
theory, that rather strong molecular weight dependences of the PS tracer diffusivities in 
PVME gels were measured, D ~ M^^ '' at gel density c = 0.200g /mL and D ~ M"^-^ at 
c = 0.235g/mL respectively^^. It was stated that the data agree well with diffusion of PS in 
PVME solution measurements where larger or equal (in the limit of large matrix molecular 
weights) diffusivities (also showing anomalous exponents) were found Clearly this aspect 



of the comparison of tracer motion in solutions and gels is explained by the PMC constraint 
release mechanism, which is not present in gels but affects the solution data depending on 
matrix (and tracer) molecular weight. Important for this experimental comparison is the 
achievement of similar matrix structural correlations in the cases of PVME solutions and 
gels^^. In PMC theory, the constraint porosity effect is important in both solution and gel 
situations and can explain the finite size fluctuation corrections. 

Figure shows PMC model calculations appropriate for the experimentally studied 
system. As argued above, in order to mimic the gel strangulation regime the specified values 



of ^g/h are chosen. From the results in figure ^ it can be seen that small variations with 
S,g/h for the other parameters are obtained. From the stated near identity of the PVME gel 
and solution density structures''^, a Flory excluded volume fractal dimension. Dp = 5/3, is 
deduced for the gels. Screening of the excluded volume interactions of the tracer is assumed 
for distances larger than the gel pore size; i.e. v = 1/2 is chosen. This choice allows us 
to connect the gel entanglement strength factor, A^, to the prefactor of the center-of-mass 
friction function, A^i; i.e. = Xo/Ne follows in eqn. (1.79). In order to crudely account for 
the known specific charge transfer attractive interactions between PS and PVME we choose 
a priori a value A^A^^e = 200, which is a factor of about 10 times larger than the value of 
Xd found in PS solutions in section 3.F. Changes of the strength factor A^ by a factor of 
2 would not result in appreciable changes of the results in figure O. From the results for 



polymer tracer diffusion in solutions shown in figure 7 of paper I, one can estimate that for 
z/ = 1/2, solution-like strength factors could lead to apparent D versus N scaling exponents 
up to 2.6. 

Most sensitive for the comparison of the PMC constraint porosity corrections to the 
measured PS in PVME tracer diffusivities is the ratio of tracer size to gel pore size spanned 
in the experiments. Following the statement in Ref. 75 we assume the characteristic pore 
length of the gel at density c = 0.235g /rriL to be = 20 A, which agrees with the size of a PS 
tracer in good solvent for a PS molecular weight of M = 5000. Note that the experiments 
span the tracer molecular weight range 10^ < M < 10^, and are therefore in the region 
where R^^^^'^^ ^ ^g. For the lower gel density, c = 0.200g /mL, we estimate the PS molecular 
weight which leads to -R^^(M) = ^p^me assuming excluded volume correlations of the 



gel, ~ c , and of the PS tracer inside the gel pores, Rg ~ M^/s for Rg < ^g] tYiis leads 
to ^g{c = 0.200g/mL) corresponding to M^^ = 5000(^)-^/^ ^ 6120. 

Figure |12| shows the PMC tracer diffusion coefficients determined with these parameters 
and thus very roughly corresponding to the experimentally measured PS diffusion coefficients 
in PVME gels^^. Whereas changes of Xg by a factor 2 and C,g/b in the employed range do not 
change the results appreciably, a larger fractal exponent u of the tracer chain, e.g. = 3/5 if 
excluded volume for the tracer was not screened, would decrease the slope of the PMC curves 
in figure |T2| down to values around D ~ A^~^'^. Changes in the fractal dimension of the gel 
only weakly affect the shape of the PMC curves. For the chosen parameters, much steeper 
molecular weight dependences than the asymptotic, reptation like law, result, and arise from 
the spatial correlations of the gel constraints. The tracer polymer feels the full entanglement 
constraints only if its size considerably exceeds the gel pore size and the length scale, b, of 
the elastic entanglement mesh. For comparison, a power law D ~ N^"^'^ as reported in Ref. 



75 is shown in figure 0. The PMC results exhibit somewhat shallower slopes. In evaluating 
the agreement with the measured tracer diffusivities one has to consider the reasonable, but 
ad hoc estimates of the parameters, the experimental uncertainties, and the complexity of 
the ternary polymer-solvent-gel system. 



5. Summary and Conclusions 

The picture of the dynamics of entangled polymeric systems as provided by the PMC theory 
shall be summarized in this final section. A comprehensive overview of all aspects of the 
PMC approach has been recently given''^. The comparisons with experimental data, dis- 
cussed in the present paper, show that an unified qualitative and sometimes quantitative, 
description of the transport properties of entangled polymer chains has been achieved. The 
comprehensive explanation of a wide variety of experimental observations is noteworthy. Al- 
ternative explanations with phenomenological theories, which exist for some of the observed 
phenomena, fail to consistently explain all of the studied experimental facts. The recent 
measurements of the dielectric relaxation times of tracer polymers in polymer melts^'^° ap- 
pear to rule out theoretical models which do not consider the collective dynamics of the 
polymer matrix. Explanations of the anomalous viscosity versus molecular weight scaling, 



rj ~ ^3.5±o.2^ provided with models considering a single polymer in a frozen environment, 
therefore, appear accidental. The same, but frequency resolved, dielectric measurements^'^° 
stress that the rapid initial decay of the disentanglement process (or its shallow high fre- 
quency loss modulus) persists even in the limit where reptation-like scaling of the final 
relaxation time is observed. This demonstrates the existence of a non-reptative effect, unre- 
lated to constraint release, which determines the primary disentanglement step. Except for 
PMC theory, where the tracer shape fluctuations are found responsible, no other consistent 
theoretical explanation apparently exists for these anomalously shallow high frequency wings 
of the dielectric and shear moduli spectra. The measurements of tracer diffusion constants 
in immobile polymer solutions, i.e. in matrices of high molecular weight polymers^^"^^, and 
in gels^^'^^, indicate that strong deviations from reptation-like asymptotes can exist even if 
constraint release corrections are absent. The observation of a ~ M~^'^ scaling in semidi- 
lute and concentrated solutions^^^^^, and of D ~ in melts^^'^^, appears to require a 
proper theoretical treatment of the polymer liquid structure. 

The quahtative picture of the polymer dynamics envisaged by PMC theory is as follows. 
Entanglements arise from the competition of polymer backbone connectivity and uncross- 
ability due to intcrmolecular excluded volume forces. Reptation-like asymptotic transport 
coefficients are obtained for polymers in fixed dense environments. This results from the 
decay of the entanglement friction via the dynamics of the collective single chain structure 
factor. Two decay channels, a coherent (center-of-mass) one for small wavevectors, and 
an internal one, termed tracer shape fluctuations, for intermediate wavevectors, are asso- 
ciated with the tracer dynamics in the friction functions. The tracer shape fluctuations 
mainly (greatly) speed up the initial disentanglement process compared to pure reptation. 
As entanglement constraints are calculated from the microscopic intcrmolecular forces, one 
material dependent strength factor, 1/a, enters into the asymptotic PMC predictions. In 
contrast to the reptation/ tube approach where a would be an universal number, in PMC 
theory 1/a measures the average, mean squared intcrmolecular force per unit density and 
therefore is non-universal in principle. Consequently, non-universal curves for the relative 
increase of the transport coefficients, e.g. 77/77^, versus rescaled molecular weight, n — N/N^., 
are predicted. Prom comparison with experiments a rather small spread in the values of a 



for some polymer melts and solutions is found, a — A ±2, and the average value lies in the 
range of the ab initio PRISM estimate. 

Finite size corrections to the transport coefficients result in PMC theory from both the 
decay of the matrix constraints, termed self-consistent constraint release, and from the 
spatial variations of their amplitudes, called constraint porosity. The spatial dependence of 
the amplitudes arises from the elastic entanglement mesh, characterized by an entanglement 
length, b, and the local compressibility correlations, specified by the density screening length 
or physical mesh, ^p. Note that in PMC theory these effects result from straightforward 
assumptions about the static structure of the polymeric liquid and its time dependence (see 
section 3.B of paper I). It can be considered a major achievement of the PMC approach that 
the variety of non-Rouse and non-reptation-like behaviors of entangled polymer dynamics 
is derived from one set of microscopically motivated approximations. This approach is 
very different from extensions of the reptation/ tube model. There new types of polymer 
chain motions are invoked in an empirical and investigator specific fashion, with parameters 
sometimes of unclear or vague relation to experimental observables, in order to account for 
unexpected (apparent) exponents characterizing the dependence of a transport coefficient on 
macromolccular size. This advantage of PMC theory is accompanied by the corresponding 
disadvantage that no intuitive predictions about the actual monomer trajectories naturally 
emerges. 

An important finding of PMC theory relates to the spatial correlations of the entan- 
glement friction in the conformational and the center-of-mass motion^^. In conceptual 
agreement with the basic reptation/ tube idea, the center-of-mass motion is determined 
from rather local friction contributions. The internal dynamics feels constraints which are 
spatially correlated across the tracer size. A first immediate consequence is the stronger 
slowing down of the internal dynamics as compared to the translational diffusion^^. An- 
other consequence are the very different effects the distinct finite size corrections have on 
the dynamics. The self-consistent constraint release mechanism, which is a more effective 
entanglement force relaxation mechanism than the tracer motion for dynamics on the Rg and 
larger scale, is important for friction contributions on global length scales. The constraint 
porosity influences the more local friction contributions most strongly. Therefore, the flrst 



finite size effect, constraint release, mainly affects the conformational dynamics, i.e. the 
internal disentanglement time and the viscosity. The constraint porosity, on the other hand, 
accelerates the center-of-mass diffusion and the end-to-end vector fluctuations of very long 
tracers. 

The constraint release mechanism is parameterized by the inverse strength factor, a, 
and therefore is rather independent of concentration and solvent quality, as follows from the 
PRISM estimates. Consequently, PMC theory explains the universal observation of apparent 
viscosity versus molecular weight power laws, f] ~ M^ "^, for both entangled polymer solutions 
and melts. The observation of a reptation-like scaling of the dielectric relaxation time in 
matrices of high molecular weight polymers supports the identification of the PMC constraint 
release mechanism as the origin of the stronger increase of the shear viscosity. Thus, the 
failure to observe the asymptotic law is predicted to be due to the lack of a clear time scale 
separation between single chain conformational relaxation and the decay of the entanglement 
forces. Such a non-Markovian situation has been invoked in the phenomenological "coupling 
model" picture of Ngai and coworkers^^'^°, and is imphcit to the "dynamic cluster" ideas of 
Douglas and Hubbard^^. 

Spatial correlations of the entanglement constraints do not play an important role in 
polymer melts. As stated and shown in figures 5 and 6 of paper I, this holds in regard to the 
internal dynamics for any polymer concentration. If translational diffusion coefficients for 
fixed reduced degrees of polymerization are considered, the small density screening length 
in melts again leads to negligible constraint porosity effects. In solutions, however, the 
center-of-mass mobility is enhanced as the amplitude of the entanglement constraints are 
proportional to the wavevector dependent osmotic compressibility. At fixed n = N/N^, the 
density screening length is larger relative to the entanglement length in solutions than in 
melts. Smaller effective friction coefficients originate from the reduced friction contributions 
on local length scales. 

Crude, a priori PRISM estimates for 6 = ^p/b lie within a factor of ^ 5 from the 
experimentally found values. In solutions, the magnitude of the uniform PMC friction co- 
efficient extracted from fits, i.e. Aij in eqns. (1.52) and (1.68), deviates from the prediction. 
Ad = 32/3a eqn. (1.53). This deviation is not understood at present, but possibly re- 



suits from our simplification of a two-component polymer-solvent system to an effective, 
one-component polymeric system. Another possible origin is the fundamental PMC ap- 
proximation, eqn. (1.24), which identifies the decay channels of the entanglement friction 
forces. The approximation to determine the decay of the entanglement constraints from the 
collective tracer structure factor, could be qualitatively correct in solutions but not quan- 
titatively accurate. Note that this approximation is supported by our findings for polymer 
melts. There, a consistent set of parameters describes the translational diffusion coefficient 
and the viscosity, as shown in section 3.C of this paper. 

Another physical situation which indicates the importance of the structural correlations 
of the matrix entanglements is the polymer tracer motion through a gel. It appears as 
support of the PMC description of the constraint porosity effects that this common mech- 
anism explains the close connection of tracer diffusivities measured in entangled solutions 
and gels^^. In our approach, (soft) gels mainly differ from entangled polymeric matrices 
in solutions because of the absence of constraint release contributions. Stronger constraint 
porosity effects than in solutions may arise depending on the fractal dimensionalities of the 
tracer and the gel, and depending on specific tracer-gel interactions. Whereas the first ef- 
fects can again be rationalized with simple considerations of the ratio of tracer size to gel 
pore size, quantitative predictions of the second effects require further experimental and 
theoretical work on polymer-gel interactions. Strong molecular weight dependences of the 
tracer diffusivities, e.g. D ~ M~^-^ in solutions^^"^^, and D ~ M~^-^ in gels^^, result from 
the enhancement of the center-of-mass motion of tracer polymers which do not sample the 
surrounding entanglement constraints uniformly as their sizes are not yet asymptotically 
large. 

Finally, we conclude with suggestions for new experiments and/ or computer simulations 
to further test out theoretical ideas about entanglements in linear chain melts, solutions, and 
gels. The key idea is to manipulate (in a controlled manner) the three fundamental length 
scales: polymer density (concentration) screening length, ^p, entanglement mesh length, b, 
and tracer radius-of-gyration, Rg. These nonuniversal structural parameters have a clear 
physical meaning, are directly measurable, and depend on chemical structure, polymer con- 
centration, temperature, solvent quality, and external pressure. Significant progress in the 



ability to test our theory against existing and future data, and to make ab initio theo- 
retical predictions, would become possible if there were more direct measurements of the 
above quantities. This is especially true in theta solvents and concentrated solutions (up 
to the melt) of not only polystyrene, but also many other materials. For the problem of 
tracer motion in gels swollen by solvent, direct measurements of the gel collective structure 
factor, tracer radius of gyration, and some reliable knowledge of the effective polymer-gel 
interaction are desirable. 

Specific fruitful areas for more experiments include the following, (i) Self and tracer dif- 
fusion (especially in the effectively frozen matrix limit) of polystyrene in good polystyrene 
solutions in the 40 — 100 % polymer concentration regime. Mapping out the concentra- 
tion and tracer degree of polymerization dependences should allow the observation of the 
crossover from reptation-like scaling of D in the melt to the "anomalous" solution scaling. 
Our theory predicts such a crossover, and can be quantitatively applied (treating the prefac- 
tor Ad, which is not understood in solutions, as a free parameter) if the required structural 
parameters are accurately known, (ii) Same as point (i), but in theta solvents covering the 
full semidilute-concentrated-melt range. Since the length scale ratio parameter 5 — ^p/h is 
larger in theta solvents (and more density dependent), we predict enhanced porosity correc- 
tions and larger effective iV-scaling exponents for the diffusion constant. {Hi) Points (i) and 
{ii) for different chemical systems. This is important because even at fixed polymer con- 
centration and solvent quality, the parameter 5 — ^p/b is nonuniversal^^ and can potentially 
be significantly varied, {iv) More tracer measurements in model gels where the probe and 
matrix (gel) polymers are identical. Alternatively, controlled variations of tracer, solvent, 
and gel chemical structure could allow the tendency for tracer adsorption or dewetting from 
the gel structure to be tuned, which in turn is predicted to affect diffusion significantly and 
in a rationally predictable manner via the parameter A^,. (v) Combined self diffusion and 
shear viscosity measurements on variable chemical systems in both the melt and solutions. 
As usual, the availability of such complementary data allows the monomer friction constant 
question to be avoided by studying the product Drj. More importantly, it places additional 
constraints on our theory by requiring the same set of parameters simultaneously explain 
the center-of-mass diffusion and conformational (stress) relaxation aspects of entangled dy- 



namics. (vi) Finally, we have discussed elsewhere^*^ our suggestions for new experiments at 
nonzero frequencies, and for transport coefficients, of entangled fluids composed of macro- 
molecules of variable fractal dimensionality (e.g. rings) and in variable spatial dimension. 
Computer simulation could potentially play a critical role here, although proper treatment 
of heavily entangled systems remains in the future. 
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Table: Material Parameters of the Different Polymer Systems Studied 



Polym. 


PBD 


PBDh 


PDMS 


PEO 


PI 


PS 


PS 


PS 


PS 


melt 
sol. 
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13 wt % 


18 wt % 


40 wt % 


Me 


1.85 10^ 


1.24 10'^ 


1.2 10"^ 


2.2 10^ 


5.0 10=^ 


1.8 10^ 


1.3 10^ 


9.6 10^ 


4.23 10^ 




31 




162 


50 


71 


173 


1250 


923 


407 


a 


4.9", 3.3'' 


4.3^= 


2.8'' 


2.6'' 


2.2^ 3.5'=' 


2.2'' 


2.8'' 


2.8'' 


2.8'' 


S 


0.01 


0.03 


0.0 


0.05 


0.05 


0.05 


0.4 


0.3 


0.4 


5(A) 


UA 


33.9 


78.6 


37.5 


62.0 


76.5 








ep(A) 


2.30 


1.75 


4.06 


1.94 


3.20 


3.95 









from viscosity; '' from diffusion; from viscosity and diffusion; from r^; ^ Mq^^^ = 28 used. 



Molecular weights of entanglement, Mg, are taken from the experimental studies whose 
data are analyzed, and entanglement degrees of polymerization, N^, are calculated from 
them using monomer repeat masses^^; entanglement lengths, b, are taken from Ref. 46 
for T = 413 K. The packing Icngths^^, p, are identified with the density screening length, 
~ 1/ gm<^^, and are also included. Melt or solution case are indicated by m or by the weight 
concentration of the polymeric component. Inverse strength parameters, a, and length scale 
ratios, S, are determined by the fits shown in figures 1 — 4, 7 and 8. 
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ure Captions 

Figure 1: Reduced viscosity, rj/rj^, for a melt of polybutadiene^^ (PBD). The solid 
line is a one-parameter PMC fit using a = 7.7, the chain curve corresponds to a = 
4.9, including the correction for the monomeric friction coefficient, Co^^^ = 2.364Co, 
suggested in Ref. 5. The dotted line is the asymptote for the second fit. The inset 
shows logiQr]/M^ shifted vertically by a constant (10) versus logarithm of reduced 
molecular weight. The lines correspond to the fits in the main figure. A vertical bar 
denotes the error bar reported in Ref. 61 arising from uncertainties in the molecular 
weight and viscosity measurements. 

Figure 2: Self diffusion constant [D, upper figure) and viscosity {t], lower figure) 
versus reduced molecular weight for a melt of hydrogenated polybutadiene^^ (PBDh) . 
The solid lines correspond to fits to D and rj with a common parameter a = 4.3; 
additionally 6 = 0.03 is found for the diffusion constant. The dotted lines are the 
asymptotic behaviors. The chain curves show the best independent fits leading to 
slightly different a (and 6 for D). The viscosity data are compared to a power law 
N^''^ (short dashes). 

Figure 3: Self diffusion constants in melts of polydimethylsiloxane (PDMS), 
polyethyleneoxide^^ (PEO), and polybutadiene^^ (PBD) versus reduced molecular 
weight. PMC fits to the individual data sets are shown using the denoted param- 
eters. For PEO and PBD, the relation = 32/3a, eqn. (1.53), is enforced. For 
PDMS, a unconstrained fit, varying a, and 6 is shown as solid line, and results in a 
small deviation of from its theoretical value, 32/3a=3.8. A power law D ~ M~^'^ 
(short dashes) is compared to the PDMS data. The high molecular weight asymptotes 
are denoted by long dashed lines. The dot-dashed curve is a fit to the PBD data with 
the best a = 4.9 from the viscosity analysis of figure 0. 

Figure 4: Self diffusion constants for polymer melts versus reduced molecular weight. 
A common fit for the data of different polymers with the indicated parameters is shown 
as a solid line. The asymptotes. Rouse and rept at ion-like, are shown with prefactors 
corresponding to the fit. 



Figure 5: Tracer diffusion coefficients versus matrix molecular weight, p — P/Pg, for 
different tracer molecular weights, n = N/N^, as specified in the figure. The model of 
Grassley^^, eqn. (1), which has been fitted to polystyrene (PS) melts by Green and 
Kramer^^ is shown by circles lor p and by squares for n — p. The PMC fit is shown 
by the solid curves; parameters are a = 3 (a priori fixed), S — 0.07 and Ad = 1.58. 
The corresponding self diffusion coefficients are indicated by a chain curve. 

Figure 6: Dielectric relaxation times, according to eqn. (2), of a polyisoprene 
tracer^ versus reduced tracer molecular weight, N/Ne, for the indicated various poly bu- 
tadiene matrix molecular weights, p = P/Pe- The lines correspond to PMC results 
from eqn. (1.70) with the parameters specified in the figure and the Table. 

Figure 7: Dielectric relaxation times for polyisoprene (PI) melts^^ versus reduced 
molecular weight. The PMC curves are calculated for the parameters a = 3.5 (solid 
line) and a = 2.2 (dashed line), where the latter value is taken from the fit to the PI 
self diffusion constants shown in figure 4 and exhibits large uncertainties; the a priori 
chosen melt-hke value 6 — 0.05 is used. 

Figure 8: Viscosities of solutions of 13 wt% (circles) and 40 wt% (diamonds) 
polystyrene in a good solvent^^'^^. The PMC fits use the parameters a — 4.0 (solid 
line) for the 13 wt% and a = 2.5 (chain curve) for the 40 wt% data. 

Figure 9: Self (open symbols) and tracer (filled symbols) diffusion constants for 
solutions of polystyrene (PS) in dibutyl phthalate (good solvent). The reported power 
law, D ~ iV~^'^, is indicated by a long dashed linc^^^^^. A common fit to the self 
(chain curve, identical tracer, N, and matrix polymer, P, degree of polymerization) 
and tracer (sohd line, P — ION) diffusion data for the three densities is shown and 
differs in the constraint release contribution only; a strong overshoot compared to the 
asymptote without finite size effects, D — D^/{1 + A^^^) (short dashes), arises from 
the constraint porosity. The theoretical connection of a and A^j is broken. 

Figure 10: Independent PMC fits to the self [S] open symbols) and tracer (T; filled 

symbols) diffusion constants for solutions of polystyrene (PS) in dibutyl phthalate 



(good solvent) from^^~^^; a — 2.8 is kept fixed for all fits, and the found parameters b 
and Ad are listed in the figure. Error bars for b — 0.33 ± 0.12, and Ad = 22 ± 10 
follow. 

• Figure 11: Low field mobilities of polystyrenesulfonate (PSS) tracers in three poly- 
acrylamide (PAA) gels of different concentration c (wt. fraction) versus tracer degree 
of polymerization, from Ref. 74. The full lines are PMC fits with the "soft gel" model, 
eqns. (1.77,1.79,1.80), with the parameters denoted; gels pore sizes, are recorded by 
stating the degree of polymerization of a PSS tracer with identical radius of gyration. 
The vertical lines mark the tracer degrees of polymerization where Rg{N) = ig{c). 

• Figure 12: Tracer diffusion constants of polystyrene (PS) as a function of tracer 
molecular weight (M) in polyvinylmcthylether (PVME) gels for parameters compara- 
ble to the study in Ref. 75 calculated from the "soft gel" model, eqns. (1.77,1.79,1.80). 
Gel pore sizes, ^g, are denoted by the molecular weights of PS tracers with identical 
radii of gyration. Ratios of pore size to entanglement length, ^/b, are used appropriate 
for polymer solutions, ~ 0.30, and hghtly strangulated gels, ^/b larger. 



A r 




_I I I I u 

12 3 4 



iog,o(N/N,) 



Fig. 1. 



12 3 




iog,o(N/N3) 

Fig. 2. 




Fig. 4. 




1 



iog,o(N/N,) 

Fig. 6. 



03 



O 



3 - 



2 - 



1 - 







PS in good solvent 




o 13 wt% 
o 40 wt % 

- a = 4.0 

- a = 2.5 







logio(N/NJ 



Fig. 8. 




logio(N/Ne) 

Fig. 9. 



o 







■2 - 




3 4 

iog,o(N) 



Fig. 11. 




V = 1 /2 Dp = 5/3 ^„N, = 200 

F g e 

= 5000 yb = 0.3 

= 5000 ijb = 0.4 

- = 5000 yb = 0.5 

- 6120 yb = 0.3 

- 6120 yb = o.4 

-^g=6120 yb = 0.5 



iog,o(M) 



Fig. 12. 



